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Abstract. We consider the non-relativistic Hartree model in the gravitational case, i.e. with 
attractive Coulomb-Newton interaction. For a given mass M > 0, we construct stationary 
states with non-zero temperature T by minimizing the corresponding free energy functional. 
It is proved that minimizers exist if and only if the temperature of the system is below a 
certain threshold T* > (possibly infinite), which itself depends on the specific choice of the 
entropy functional. We also investigate whether the corresponding minimizers are mixed or 
pure quantum states and characterize a critical temperature T c £ (0, T*) above which mixed 
states appear. 



1. Introduction 

In this paper we investigate the non-relativistic gravitational Hartree system with temper- 
ature. This model can be seen as a mean-field description of a system of self-gravitating 
quantum particles. It is used in astrophysics to describe so-called Boson stars. In the present 
work, we are particularly interested in thermal effects, i.e. (qualitative) differences to the zero 
temperature case. 

A physical state of the system will be represented by a density matrix operator p € 
©i(L 2 (R 3 )), i.e. a positive self-adjoint trace class operator acting on L 2 (M 3 ;C). Such an 
operator p can be decomposed as 

(1.1) p = J2 x o\^)(^\ 

with an associated sequence of eigenvalues (Xj)j^ G I 1 , Xj > 0, usually called occupation 
numbers, and a corresponding sequence of eigenfunction forming a complete or- 

thonormal basis of L 2 (R 3 ), cf. [33J. By evaluating the kernel p(x,y) on its diagonal, we 
obtain the corresponding particle density 

n p (x) = ^A J |^(x)| 2 G4(M 3 ). 

In the following we shall assume that 

(1.2) / n p (x)dx = M, 
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for a given total mass M > 0. We assume that the particles interact solely via gravitational 
forces. The corresponding Hartree energy of the system is then given by 

£h[p] :=£kin[p] -£ P otM =tr(-Ap) --tr(V p p) , 
where V p denotes the self-consistent potential 

V p = n p * — 

and V is the usual convolution w.r.t. x S M 3 . Using the decomposition (jl.ip for p, the Hartree 
energy can be rewritten as 

%[/>]=!>/ |V^(x)| 2 dx-i ff n rW n »M dxdy 

To take into account thermal effects, we consider the associated free energy functional 
(1.3) T T [p}:=S H [p}-TS[p] 
where T > denotes the temperature and S[p] is the entropy functional 

S[p] :=-trf3(p) . 

The entropy generating function (3 is assumed to be convex, of class C 1 and will satisfy some 
additional properties to be prescribed later on. The purpose of this paper is to investigate 
the existence of minimizers for Tt with fixed mass M > and temperature T > and study 
their qualitative properties. These minimizers, often called ground states, can be interpreted 
as stationary states for the time-dependent system 

(1-4) i ^P(t) = [H p (t) , p(t)] , p(0) = p^ • 

Here [A, B] = AB — B A denotes the usual commutator and H p is the mean-field Hamiltonian 
operator 

(1.5) H p :=-A-n p * ^- . 

Using again the decomposition (jl.ip . this can equivalently be rewritten as a system of (at 
most) countably many Schrodinger equations coupled through the mean field potential V p : 

j i d t ipj + V(t, x) ipj = , j'GN, 

\ -AV P = 4^. eN A,|^(^)| 2 - 

This system is a generalization of the gravitational Hartree equation (also known as the 
Schrodinger- Newton model, see [5]) to the case of mixed states. Notice that it reduces to a 
finite system as soon as only a finite number of Xj are non-zero. In such a case, p is a finite 
rank operator. 

Establishing the existence of stationary solutions to nonlinear Schrodinger models by means 
of variational methods is a classical idea, cf. for instance [15]. A particular advantage of such 
an approach is that in most cases one can directly deduce orbital stability of the stationary 
solution w.r.t. the dynamics of (jl.4p or, equivalently, (jl.6p . In the case of repulsive self- 
consistent interactions, describing e.g. electrons, this has been successfully carried out in 
El [U [23]. In addition, existence of stationary solutions in the repulsive case has been 
obtained in [231 E3 EH EZ] using convexity properties of the corresponding energy functional. 



(1.6) 
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In sharp contrast to the repulsive case, the gravitational Hartree system of stellar dynamics, 
does not admit a convex energy and thus a more detailed study of minimizing sequences is 
required. To this end, we first note that at zero temperature, i.e. T = 0, the free energy 
J~t[p] reduces to the gravitational Hartree energy £h[p\- For this model, existence of the 
corresponding zero temperature ground states has been studied in [14|. [171 US] and, more 
recently, in [5]. Most of these works rely on the so-called concentration- compactness method 
introduced by Lions in [18J. According to [J3], it is known that for T = the minimum of 
the Hartree energy is uniquely achieved by an appropriately normalized pure state, i.e. a rank 
one density matrix po = M |V ; o)(V'o|- The concentration-compactness method has later been 
adapted to the setting of density matrices, see for instance [13] for a recent paper written this 
framework, in which the authors study a semi-relativistic model of Hartree- Fock type at zero 
temperature. 

Remark 1.1. In the classical kinetic theory of self-gravitating systems, a variational approach 
based on the so-called Casimir Junctionals has been repeatedly used to prove existence and 
orbital stability of stationary states of relativistic and non-relativistic Vlasov-Poisson models: 
see for instance 03H35H3Sll2Sll22l32lEll3iIlEIl- These functionals can be regarded as the 
classical counterpart of Ft [p] and such an analogy between classical and quantum mechanics 
has already been used in [2J1 [7J E] • 

In view of the quoted results, the purpose of this paper can be summarized as follows: First, 
we shall prove the existence of minimizers for Ft, extending the results of [14\ \T7\ \19\ [5] to 
the case of non-zero temperature. As we shall see, a threshold in temperature arises due to the 
competition between the Hartree energy and the entropy term and we find that minimizers of 
Tt exist only below a certain maximal temperature T* > 0, which depends on the specific form 
of the entropy generating function /3. One should note that, by using the scaling properties of 
the system, the notion of a maximal temperature for a given mass M can be rephrased into 
a corresponding threshold for the mass at a given, fixed temperature T. Such a critical mass, 
however, has to be clearly distinguished from the well-known Chandrasekhar mass threshold 
in semi-relativistic models, cf. [TBJ [TTJ [TH] . Moreover, depending on the choice of /3, it could 
happen that T* = +oo, in which case minimizers of Tt would exist even if the temperature 
is taken arbitrarily large. In a second step, we shall also study the qualitative properties of 
the ground states with respect to the temperature T E [0, T*). In particular, we will prove 
that there exists a certain critical temperature T c > 0, above which minimizers correspond to 
mixed quantum states, i.e. density matrix operators with rank higher than one. If T < T c , 
minimizers are pure states, as in the zero temperature model. 

In order to make these statements mathematically precise, we introduce 



The set fj can be interpreted as the cone of nonnegative density matrix operators with finite 
energy. Using the decomposition (jl.ip . if p G $), we obtain that ipj G if 1 (M 3 ) for all j £ N 
such that Xj > 0. Taking into account the mass constraint (|1.2p we define the set of physical 
states by 




} 



and consider the norm 




£m :={p££ : trp = M} . 
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We denote the infimum of the free energy functional Ft, defined in (|1.3p . by 

(1.7) i M ,T = inf F T [p] ■ 

The set of minimizers will be denoted by SDTm C $)m- As we shall see in the next sec- 
tion, im,t < if DJIm 7^ 0- This however can be guaranteed only below a certain maximal 
temperature T* = T*(M) given by 

(1.8) T*(M) : = sup{T > : i M ,T < 0} . 

This maximal temperature T* will depend on the choice of the entropy generating function 
(3 for which we impose the following assumptions: 

(f31) (3 is strictly convex and of class C 1 on [0,oo), 
(/32) f3 > on [0, 1] and 0(0) = /3'(0) = 0, 

{(33) sup me(0jOo) < 3. 

A typical example for the function (3 reads 

P(s)= s p, pe(l,3]. 

Such a power law nonlinear ity is of common use in the classical kinetic theory of self- 
gravitating systems known as polytropic gases. One of the main features of such models is to 
give rise to orbitally stable stationary states with compact support, cf . [lOj EHl GQj EH G21 ES] , 
clearly a desirable feature when modeling stars. We shall prove in Section 6, that T* is finite 
Hp is not too large. The limiting case as p approaches 1 corresponds to /3(s) = s Ins but in 
that case the free energy functional is not bounded from below, see |21j for a discussion in 
the Coulomb repulsive case, which can easily be adapted to our setting. 

Up to now, we have made no distinction between pure states, corresponding density matrix 
operators with rank one, and mixed states, corresponding to operators with finite or infinite 
rank. In [TJ] Lieb has proved that for T = minimizers are pure states. As we shall see, this is 
also the case when T is positive but small and as a consequence we have: %m,t = i>M,o+T (3(M). 
Let us define 

(1.9) T C (M) := max {T > : i M ,T = iufi + r P[M) V r G (0, T] } . 

With these definitions in hand, we are now in the position to state our main result. 

Theorem 1.1. Let M > and assume that (/31)-(/33) hold. Then, the maximal tempera- 
ture T* defined in (j 1 . 8|) is positive, possibly infinite, and the following properties hold: 

(i) For all T < T* , there exists a density operator p £ Sjm such that Ft[p] = im,t- 
Moreover p solves the self- consistent equation 

p = {f3'r\{»-H p )/T) 

where H p is the mean-field Hamiltonian defined in (jl.5p and p < denotes the La- 
grange multiplier associated to the mass constraint. 

(ii) The set of all minimizers 9JIm C^m is orbitally stable under the dynamics of (|1.4p . 

(iii) The critical temperature T c defined in (|1.9p is finite and a minimizer p G 2Rm is a 
pure state if and only ifTG [0, T c ] . 

(iv) //, in addition, (3(s) = s p with p £ (1, 7/5), then T* < +oo. 
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The proof of this theorem will be a consequence of several more detailed results. We shall 
mostly rely on the concentration-compactness method, adapted to the framework of trace 
class operators. Our approach is therefore similar to the one of [6] and |13| . with differences 
due, respectively, to the sign of the interaction potential and to non-zero temperature effects. 
Uniqueness of minimizers (up to translations and rotations) is an open question for T > T c . 
For T G [0, T c ] , the problem is reduced to the pure state case, for which uniqueness has been 
proved in [H] (also see [T2]). 

This paper is organized as follows: In Section [2] we collect several basic properties of the 
free energy. In particular we establish the existence of a maximal temperature T* > and 
derive the self-consistent equation for p G $)m- I n Section [31 we derive an important a 
priori inequality for minimizers, the so-called binding inequality, which is henceforth used in 
proving the existence of minimizers in Section SJ Having done that, we shall prove in Section 
[5] that minimizers are mixed states for T > T c , and we shall also characterize T c in terms 
of the eigenvalue problem associated to the case T = 0. In Section [6J we shall prove that 
T* is indeed finite in the polytropic case, provided p < 7/5 and furthermore establish some 
qualitative properties of the minimizers as T — > T* < +oo. Finally, Section [7] is devoted to 
some remarks on the sign of the Lagrange multiplier associated to the mass constraint and 
related open questions. 

2. Basic properties of the free energy 

2.1. Boundedness from below and splitting property. As a preliminary step, we ob- 
serve that the functional Tt introduced in (|1.3|) is well defined and im,t > — oo. 

Lemma 2.1. Assume that (f31)-(f32) hold. The free energy Tt is well-defined on Sjm an d 
iu,T is bounded from below. If J~t[p] is finite, then sJrTp is bounded in H 1 (M. 3 ). 

Proof. In order to establish a bound from below, we shall first show that the potential energy 
£pot[/o] can be bounded in terms of the kinetic energy. To this end, note that for every p £ ft 
we have 

£ P ot[p] < CllnpH^ 2 llnpH^ 2 

by the Hardy-Littlewood-Sobolev inequality. Next, by Sobolev's embedding, we know that 
1 1 rip 1 1^3 is controlled by HVy^H^ which, using the decomposition (jl.ip . is bounded by 
tr(— A p). Hence we can conclude that 

(2.1) £ pot [p] < C\\n p \\% 2 tr(-Ap) 1 / 2 

for some generic positive constant C. By conservation of mass, the free energy is therefore 
bounded from below on Sjm according to 

T T [p] > tr(-Ap) - CM 3/2 tr(-Ap) 1/2 > -j C 2 M 3 

uniformly w.r.t. p G Hmi thus establishing a lower bound on «m,t- For the entropy term 
S[p] = — tr f3(p) we observe that, since is convex and /3(0) = 0, it holds < /?(/?) < /3(M) p 
for all p G S) and (3(p) G ©i, provided p G ©i. Hence, all quantities involved in the definition 
of Ft are well-defined and bounded on Sjm- D 

Throughout this work, we shall use smooth cut-off functions defined as follows. Let x De 
a fixed smooth function on R 3 with values in [0,1] such that, for any x G M 3 , xi x ) = 1 if 
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\x\ < 1 and x( x ) = if \x\ > 2. For any R > 0, we define \R an d £r by 



(2.2) X r{x) = x{x/R) and £ R (x) = y/l- x{x/R) 2 V x G R 3 . 

The motivation for introducing such cut-off functions is that, for any u G // 1 (IR 3 ) and any 
potential V, we have the identities 

/ |u| 2 dx = / \xRu\ 2 dx+ / \^ R u\ 2 dx 

JR3 Vr3 Jr3 

and / F|u| 2 dx= / V \xr u\ 2 dx + / V\^Ru\ 2 dx, 

JW 3 7 R 3 J R 3 

and the IMS truncation identity 

(2.3) I \V( X Ru)\ 2 dx+ [ \X7(£ RU )\ 2 dx= [ \X7u\ 2 dx- [ \u\ 2 V • (V X r + V£*) dx . 

JR3 J R 3 J R 3 J R 3 V w ' 

=0(ii- 2 ) as _R->oo 

A first application of this truncation method is given by the following splitting lemma. 
Lemma 2.2. For p G Sjm, we define p^ = XRPXR. an d Pr = £,rp£,r- Then it holds: 

S[p { £] + S[p (2) ] > S[p] and fjdn^g 5 ] + £kin[/4 2) ] < £kin[p] + 0{R- 2 ) asR^+oo. 
Proof. The assertion for ^kinl/ 5 ] is a straightforward consequence of (|2.3p . namely 

t^-Ap^) + tr(-Apg } ) = tr(-Ap) + 0(iT 2 ) as R -> +oo . 
For the entropy term, we can use the Brown-Kosaki inequality (cf. [2]) as in [6j Lemma 3.4] 



to obtain 



tr^^+tr^pg^tr^p). 



□ 



2.2. Sub-additivity and maximal temperature. In order to proceed further, we need 
to study the dependence of iyiT with respect to M and T and prove that the maximal 
temperature T* as defined in (|1.8p is in fact positive. To this end, we rely on the translation 



invariance of the model. For a given y G M 3 , denote by r y : L 2 (M 3 ) — > L 2 (M 3 ) the translation 
operator given by 

(V) = /(•-!/) V/GL 2 (M 3 ). 

Proposition 2.3. Let %m,t be given by (|1.7p and assume that (f31)-(f32) hold. Then the 
following properties hold: 

(i) As a function of M , imt is non-positive and sub-additive: for any M > 0, m G (0, M) 
and T > 0, we have 

lM,T < Ul-m,T + im,T < . 

(ii) The function %m,t is a non-increasing function of M and a non- decreasing function 
of T . For any T > 0, we have im,t < if and only ifT<T*. 

(hi) For any M > 0, T*(M) > is positive, possibly infinite. As a function of M it is 
increasing and satisfies 

fry 

T*(M) > max —— \i l)0 \ . 
0<m<M j3(m) 

As a consequence, T* > and T*(M) = +oo for any M > i/lim< j _ > o + P(s)/s 3 = 0. 
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Proof. We start with the proof of the sub-additivity inequality. Consider two states p G SjM-m 
and a G S) m , such that Jr[p] < iM-m,T + e and J-t[c] < i m< T + £■ By density of finite rank 
operators in 35 and of smooth compactly supported functions in L 2 , we can assume that 

J 

p = ^Z x j\^j)^j\ » 

with smooth eigenfunctions (ipj)j =1 having compact support in a ball -B(0, R) C R 3 , for some 
J G N. After approximating <r analogously, we define cr# e := T^ Re arsR e , where e 6 § 2 C K 3 
is a fixed unit vector and r is the translation operator defined above. Note that we have 
P o~R e = 0\R e p = 0, hence p + <7Re G i^M and tr /3(p + o-R e ) = tr f}(p) + tr ft (a Re)- Thus we have 

*M,T < ^r[p + O-Re] = F T [p] + Ft[<A + 0(1/ R) < lM-m,T + i m ,T + 2e , 

where the 0(1/R) term has in fact negative sign so that we can simply drop it. Taking the 
limit e — > yields the desired inequality. 

Next, consider a minimizer p of 8r subject to tr p = M. It is given by an appropriate 
rescaling of the pure state obtained in [T3]. For an arbitrary A G (0,oo), let (U\f)(x) := 
\ 3 / 2 f(\x) and observe that p\ := U^pU\ G Sjm- As a function of A, the Hartree energy 
£h[px] = A 2 <Ski n [p] — A Spot [p] has a minimum for some A > 0. Computing ^£h[p\] = 0, we 
infer that A = £ po t[p]/(2£kin[p]) and moreover 

. l (gpotM) 2 

As a consequence, we have zm,o = ^1,0 and 

(2.4) F T [p} = i M ,o + Tp(M) = f3(M) (t - ^yl*i,o|) > *m,t , 

thus proving that ^m,t < for T small enough. 

Since f3 is non-negative function on [0, 00), the map T 1— > J~t[p\ is increasing. By taking 
the infimum over all admissible p G SjMi we infer that T 1— > im,t is non-decreasing. The 
function M \-t im,t is non-increasing as a consequence of the sub-additivity property. As a 
consequence, T*(M) is a non-decreasing function of M, such that 

T*(M) > lim T*(M) . 

By the sub-additivity inequality and (|2.4|) . we obtain 

M 3 / M 3 \ 

iM,T<ni M/n;T <nP(f)T--^ M=n/3(f) I T - T^M) I* 1 ' ' ) 

for any n G N*. Since lim s _>.o + P(s)/s = 0, we find that im,t < by passing to the limit as 
n — > 00. In the particular case lim s _ > .o + /5(s)/s 3 = 0, we conclude that T*(M) = +00 for any 



M > 0. Similarly, using again the sub-additivity inequality and (|2.4p . we infer 



m 



im,t < i m ,T < (3{m) T - — — |*i,o I V m G (0, M] , 

V Pi™) J 

which provides the lower bound on T*(M) in assertion (iii). By definition of T*(M), we also 
know that iu,T is negative for any T < T*(M). From the monotonicity of T h-> im,t> we obtain 
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that %m,t = if T > T* and T* < oo. Because of the estimate zm,t < ^m,t + (T — To) (3(M) 
for any T > Tq, we also find that %m,t* = if T* < oo. □ 

2.3. Euler-Lagrange equations and Lagrange multipliers. As in [HJ [6] , we obtain the 
following characterization of p G 9Hm- 

Proposition 2.4. Lei M > 0, T £ (0,T*(M)] and assume that fl31)-((32) hold. Consider a 
density matrix operator p G Sjm which minimizes Tt ■ Then p is such that 

(2.5) tr(V» = 4tr(-Ap) 
and satisfies the self- consistent equation 

(2.6) p=(PT 1 ((p-H p )/T), 

where H p is the mean-field Hamiltonian defined in (|1.5|) and p < denotes the Lagrange 
multiplier associated to the mass constraint tr p = M. Explicitly, p is given by 

(2-7) p = ±tv((H p + Tf3'(p))p) . 

Proof. Let p G 9Jtj;/ be a minimizer of J-^- Consider the decomposition given by (jl.ip . If we 
denote by the density operator in S)m given by 

PA = A 3 ^A i |^(A-))(^(A-)| , 

then, as in the proof of Proposition 12.31 we find that £h[p\] = ^ 2 £km[p] — A£ po t[/?] while 
S[p\] = S[p] for any A > 0. Hence the condition -^£h[p\\\\=i = exactly amounts to 
£pot[p] = 2£kin[p]- Next, let a G S)m- Then (1 - t) p + t a G f)Af and 

t h-> J" T [(1 -i)/9 + icr] 

has a minimum at t = 0. Computing its derivative at t = and arguing by contradiction 
implies that p also solves the linearized problem 

inf tr {(H p + Tf3'(p))(a-p)) . 

Computing the corresponding Euler-Lagrange equations shows that the minimizer of this 
problem is p = (/3')~ 1 (( / u — H p )/T) where p denotes the Lagrange multiplier associated to 
the constraint tr p = M. Since the essential spectrum of H p is [0, oo), we also get that p, < 
since p is trace class and (/3') _1 > on (0, oo). □ 

Using the decomposition we can rewrite the stationary Hartree model in terms of (at 
most) countably many eigenvalue problems coupled through a nonlinear Poisson equation 

J Ai/jj + V p iPj = , j G N , 

\ -AU p = 47r^. eN A i |V i | 2 , 

where (/ij)j6N £ K denotes the sequence of the eigenvalues of and (ipj,ipk)i 2 = ^',fc- 
The self-consistent equation (j2.6|) consequently implies the following relation between the 
occupation numbers (Aj)jgN and the eigenvalues (pj)j £ ^: 

(2.8) A i = 03 , )- 1 (( M -^)/T) +J 

where s + = (s + |s|)/2 denotes the positive part of s. Upon reverting the relation (|2.8p we 
obtain pj = p — T f3'{\j) for any pj < p. 
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The Lagrange multiplier p is usually referred to as the chemical 'potential. In the existence 
proof given below, it will be essential, that p < 0. In order to show that this is indeed the 
case, let p(M) := sup m6(0)M] ffi . If p G Sj M , then 

tr(/3'(p)p) <p(M) tr/3(p) . 
Notice that if (/33) holds, then p(M) < 3. 

Lemma 2.5. Let M > and T < T*(M). Assume that p E Sjm is a minimizer of Tt and 
let p be the corresponding Lagrange multiplier. With the above notations, if p{M) < 3, then 
M n<p(M)i M ,T < 0. 



Proof. By definition of %m,t and according to (|2.7p , we know that 

i M ,T = ^{-/\p-\V pP + Tj3{p)) , 

Mp = tv(-Ap-V pP + Tp'(p)p) . 
Using (|2.5p , we end up with the identity 

p(M) i M ,T - Mix = (3 - p(M)) tr(-A p) + T tr (p(M) /3(p) - /3'(p) p) > , 
which concludes the proof. □ 

The negativity of the Lagrange multiplier p, is straightforward in the zero temperature 
case. In our situation it holds under Assumption (/33), but has not been established for 
instance for f3(s) = s p with p > 3. In fact, it might even be false in some cases, see Section [7] 
for more details. 

Corollary 2.6. Let T > 0. Then M i-> iu,T is monotone decreasing as long as T < T*(M) 
and p(M) < 3. 

Proof. Let p E Sjm be such that 7"t[/o] < im,t + for some e > to be chosen. With no 
restriction, we can assume that £ po t[p] = 2<?kin[/ 3 ] and define p[p] := ^ .Ft [A p]|a=i- The same 
computation as in the proof of Lemma 12.51 shows that 

p(M) (i M ,T + e) - Mp, > (3 - p(M)) tr(- A p) + T tr (p(M) - p'{p) p) > , 

since, by assumption, p(M) < 3. This proves that M/i[p] < im,t/2 < for any e G 
(0, |iA/,T|/2), if p(M) < 3. This bound being uniform with respect to p, monotonicity easily 
follows. □ 



Remark 2.7. Under the assumptions of Lemma |2.5| we observe that 

d 

d\ 

provided p(M) < 3 and p £ $)m, which proves the strict monotonicity of M \- > im,t- However, 
at this stage, the existence of a minimizer is not granted and we thus had to argue differently. 



T T [\ P ]\x=i = pM <0, 



3. The binding inequality 

In this section we shall strengthen the result of Proposition 12.31 (i) and infer a strict sub- 
additivity property of iM,T, which is usually called the binding inequality; see e.g. |13j . This 
will appear as a consequence of the following a priori estimate for the spatial density of the 
minimizers. 
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Proposition 3.1. Let p S Sjm be a minimizer of J~t- There exists a positive constant C such 
that, for all R > sufficiently large, 

[ n p (x) dx < ^ . 
J\x\>R K 

This result is the analog of |13l Lemma 5.2]. For completeness, we shall give the details of 
the proof, which requires p < 0, in the appendix. The following elementary estimate will be 
useful in the sequel. 

Lemma 3.2. There exists a positive constant C such that, for any p £ S)m> 



f !hM d x<CM 3 / 2 (tr(-Ap)) 1 



Proof. Up to a translation, we have to estimate J* K3 n p (x) dx and it is convenient to split 
the integral into two integrals corresponding to \x\ < R and |x| > R. By Holder's inequality, 
we know that, for any p > 3/2, 

n (x) , / „_i \(p-i)/p „ _ 2 P -3 



'Br l x l x ' 

where Br denotes the centered ball of radius R. Similarly, for any p < 3/2, 

n (x) , / _i \(p-i)/p „ , _2p-3 
c ^d*<(4^) ||n p || iPJ R ^ . 

Applying these two estimates with, for instance, p = 3 and p = 6/5 and optimizing w.r.t. R > 
0, we obtain a limiting case for the Hardy-Littlewood-Sobolev inequalities after using again 
Holder's inequality to estimate ||n p || L 6/5 in terms of ||n p ||£i and ||n p ||^3: 

n p( x ) ||3/2,, |,l/2 



. p\\ L i \\>"p\\ L z ■ 

R 3 l x l 

We conclude as in (|2.ip using Sobolev's inequality to control ||n p ||/,3 by tr(— Ap). □ 

As a consequence of Proposition 13.11 and Lemma 13.21 we obtain the following result. 

Corollary 3.3 (Binding inequality). Let jVfW > and > 0. If there are minimizers 

for ij^(i) T and ijy(a) j>, i/ien 

Proof. Consider two minimizers and p( 2 ) for T and « M (2) T respectively and let xr 
be the cut-off function given in (|2.2p . By Lemma 12.21 we have 

tr(-A ( W> < tr(-A /0 W) + O(i?- 2 ) and tr/3( XR pW Xfl ) < tr/3(p«) . 
To handle the potential energies, we observe that 



£ P ot[XRP {e) XR]-£ P ot[p (l) } < II ' <^~>~y°» -p EZ^Zl dxdy 



i - x 2 r( x ) x 2 R (y)) n P w 0) n P w (y) 



x-y\ 



n p {t){x)n p{ i) (y) 



<l , | dxdy . 

l{\x\>R}x{\y\>R} \X-y\ 
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Using Lemma 13.11 and Lemma 13.21 we obtain 

r -i 1/2 r 

£MXRP {£) XR] - £ P ot[p {i) ) < C tr(-A/)) ' / n w {x) dx < 0(R- 2 ) 

1 J J\x\>R 

for R > large enough. This shows that, for any R > sufficiently large 

FtIxr P {t> XR] < i M W,T + 0(R- 2 ) for £ = 1,2. 
Consider now the test state 

PR-=XR XR + T^Re XR P (2) XR T 5 Re 

for some unit vector e 6 S 2 . Since Hn^H^i < + M^ 2 \ by monotonicity of M \— > %m,t 
(see Proposition 12.31 (ii)). we get 

(1\ (2) MMMW 
i M (i)+M(2),T ^ ^t[pk] < Ft[xrP k ' Xr] +J 7 t[xrP { ' Xr] g-^ 

^ Z M(!),T + Z M( 2 ),T + ^2 9^ 

for some positive constant C, which yields the desired result for R sufficiently large. □ 

4. Existence of minimizers below T* 

By a classical result, see e.g. [1'6\ Corollary 4.1], conservation of mass along a weakly 
convergent minimizing sequence implies that the sequence strongly converges. More precisely, 
we have the following statement. 

Lemma 4.1. Let (pk)ken £ £)m be a minimizing sequence for Ft, such that pt — p weak—* 
in and n Pk — > n p almost everywhere as k — )■ oo. TTien Pk ^ P strongly in S) if and only if 
trp = M. 

Proof. The proof relies on a characterization of the compactness due to Brezis and Lieb (see [1] 
and |15} Theorem 1.9]) from which it follows that 



lim ( / n Pk dx — \n p — n Pk \ dx \ = n p dx 



k 

and lim ( tr(-Ap) - tr ( - A (p - p k )) ) = tr(-Ap) . 

k— 5>oo V / 

By semi-continuity of Tt% monotonicity of M \— > im,t according to Proposition 12.31 (ii) and 
compactness of the quadratic term in £h, we conclude that lim^oo tr(— A (p — pf.)) = if 
and only if tr p = M. □ 

With the results of Section [2] in hand, we can now state an existence result for minimizers 
of Ft- To this end, consider a minimizing sequence (p n )n<=N for Tt and recall that (p n )neN is 
said to be relatively compact up to translations if there is a sequence (a n ) ng N of points in M 3 
such that r* n p n T an strongly converges as n — > oo, up to the extraction of subsequences. 

Clearly, the sub-additivity inequality given in Lemma 12.31 (i) is not sufficient to prove the 
compactness up to translations for (p n )n&n- More precisely, if equality holds, then, as in the 
proof of Lemma [2. 3\ one can construct a minimizing sequence that is not relatively compact in 
S) up to translations. This obstruction is usually referred to as dichotomy, cf. [18]. To overcome 
this difficulty, we shall rely on the strict sub-additivity of Corollary 13.31 which, however, only 
holds for minimizers. This is the main difference with previous works on Hartree-Fock models. 



12 G.L. AKI, J. DOLBEAULT, AND C. SPARBER 

As we shall see, the main issue will therefore be to prove the convergence of two subsequences 
towards minimizers of mass smaller than M. 

Proposition 4.2. Assume that (J31)-((38) hold. Let M > and consider T* = T*{M) defined 
by (|1.8|b For all T <T* , there exists an operator p in Sjm such that J^tIp] = ihi,T- Moreover, 
every minimizing sequence (p n ) n m for iu,T is relatively compact in S) up to translations. 



Proof. The proof is based on the concentration-compactness method as in [13] . Compared to 
previous results (see for instance [201 ED 1221 02] ) , the main difficulty arises in the splitting 
case, as we shall see below. 

Step 1: Non-vanishing. We split 

£M=H n -^M dxdy 
JJrs \x-y\ 

into three integrals I±, I2 and ^3 corresponding respectively to the domains \x — y\ < 1/R, 
1/R < \x — y\ < R and \x — y\ > R, for some R > 1 to be fixed later. Since n Pn is bounded 
in L 1 (M 3 ) n L 3 C L 7 / 5 (R 3 ) by Lemma \2.1\ by Young's inequality we can estimate I\ by 

.. n2 hi 1 \ || C 

h < \\n Pn \\ L r/5 III • I WlV\b 1/r ) < ^77 » 
and directly get bounds on I2 and I3 by computing 



h<R n Pn (x) n Pn (y) dxdy <RM sup / 

JJ\x-y\<R y eM?Jy+B R 



n Pn {x) dx 



1 f f M 2 

h < J J n Pn (x) n Pn (y) dxdy < — . 

Keeping in mind that i m,t < 0, we have 

FtIpti] > ihi.T > -h - h - h 
for any n large enough, which proves the non-vanishing property: 

f / \ i 1 ( M 2 C \ 

-~ Lb r npJx) dx ~RM{- W ~ -R ~ l^J > ° 

for R big enough and for some sequence (a n ) n &N of points in M 3 . Replacing p n by r* n p n T an 
and denoting by p^ the weak limit of (p n )neN (up to the extraction of a subsequence), we 
have proved that M ^ = J" R3 n (1) dx > 0. 

Step 2: Dichotomy. Either M^' = M and p n strongly converges to p in f) by Lemma 14. 1| 
or M"W G (0,M). Let us choose R n such that J R3 n (1) dx = + (M — M^)/n where 

Pn ■= XRn Pn XR„- Let p^ := ^ /3 n By definition of R n , Hindoo R n = 00. By Step 1, 
we know that pn^ strongly converges to p^K By Identity (|2.3p and Lemma 12.21 we find that 

r [■ n (i) (x) n ( 2 ) (y) 
F T [ Pn \>F T [p^]+F T [p^] + 0(R- 2 )- // ^ dxdy, 

JJr 3 xr 3 \ x -y\ 

thus showing that 

iM,T = hm T T [Pn] > ^r[p (1) ] + hm Jr[/4 2) ] . 
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By step 1, linin^oo J* R3 n (2) dx = M — M^> . By sub-additivity, according to Lemma I2T31 (i), 
p^ is a minimizer for i M (i) T , ( P n )neN is a minimizing sequence for i M _ M (i) T and 

«M,T = HlW,T + Hl-MW,T ■ 

Either i M _ M (i) T = and then zm,t = ^M-itfW t» which contradicts Corollary 12. 6} and the 
assumption T < T*, or i M _ M (i) T < 0. In this case, we can reapply the previous analysis 

to (pn jneN and get that for some > 0, (/9n )«eN converges up to a translation to a 

minimizer p^ 2 ** for ijvf( 2 ) r an d 

«M,T = «M(!),T + *A/( 2 ),T + ^M-MW-M( 2 ),T • 

From Corollary 13.31 and 12.31 (i), we get respectively g'mW+mW t < ^Mt 1 ) r + ^m( 2 ) t an< ^ 
«m«+m( 2 ),t + *m-m(i)-m( 2 ),t ^ »Af,T, a contradiction. □ 

As a direct consequence of the variational approach, the set of minimizers 9Hm is orbitally 
stable under the dynamics of (|1.4p . To quantify this stability, define 

dista^cr) := inf ||p-<r||jj . 

Corollary 4.3. Assume that fl31)-(P3) hold. For any given M > 0, let T G (0,T*(M)). For 
any e > 0, i/iene exists <5 > suc/i that, for all p in £ F)m with distoj{ M (pi n ) < 8, 

sup dist mM (p(t)) < e 

where p(t) is the solution of (|1.4p with initial data p m G i^Af- 

Similar results have been established in many earlier papers like, for instance in [23J in 
the case of repulsive Coulomb interactions. As in [U [23], the result is a direct consequence 
of the conservation of the free energy along the flow and the compactness of all minimizing 
sequences. According to [H], for T G (0, T c ], the minimizer corresponding to im,t is unique 
up to translations (see next Section). A much stronger stability result can easily be achieved. 
Details are left to the reader. 

5. Critical Temperature for mixed states 

In this subsection, we shall deduce the existence a critical temperature T c G (0,T*), above 
which minimizers p G 9Ha/ become true mixed states, i.e. density matrix operators with rank 
higher than one. 

Lemma 5.1. For all M > 0, the map T 1— >• im,T * s concave. 
Proof. Fix some Tq > and write 

T T [p}=J r T [p] + (T-T )\S[ P ]\ . 
Denoting by P t the minimizer for J-t , we obtain 

iM,T < iM,T + {T- T ) \S[pr ]\ 

which means that |<S[pt ]| lies in the cone tangent to T 4 im,t and iu,T lies below it, 
i.e. T 1— > im T is concave. □ 
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Consider T c defined by (|1.9p . i.e. the largest possible T c such that im,t = im,o + T f3(M) 
for T E [0, T c ] and recall some results concerning the zero temperature case. Lieb in [T3] 
proved that Tt=o = £h has a unique radial minimizer po = M \ipo)(ipo\- The corresponding 
Hamiltonian operator 



(5.1) H :=-A-\i> \ 2 *\-\- l = H f 

admits countably many negative eigenvalues (p^)j^, which accumulate at zero. We shall 
use these eigenvalues to characterize the critical temperature T c . To this end we need the 
following lemma. 

Lemma 5.2. Assume that (/31)-(/33) hold. With T c defined by (fL9j) . T C (M) is positive for 
any M > 0. 

Proof. Consider a sequence (T n ) ng N £ such that lim^oo T n = 0. Let p( n ) E 9#Af denote 
the associated sequence of minimizers with occupation numbers (Xj)j^. According to (|2.8|) . 
we know that 

AfM^r 1 ((/>M n VTn) VjeN, 

where, for any n E N, (^^ )j 6 n denotes the sequence of eigenvalues of JET^n) and p^ < 
is the associated chemical potential. Since p( n ' is a minimizing sequence for J-t=q, we know 
that 

lim /4 n) = /x? < 

where (p,j)j & ^ are the eigenvalues of i^o- Arguing by contradiction, we assume that 

liminf Al n) = e > . 

n— >oo 

By (|2.8p and the fact that f3' is increasing, this implies: p^ > p(^ — > p® as rt — > oo. Then 
M = An > lim A^ n) = lim^')- 1 ( ^^-) > lim^T 1 ( f M^ > ) = +oo . 

This proves that there exists an interval [0, T c ) with T c > such that, for any T n E [0,T C ), it 
holds // n ) < p^, and, as a consequence, p^ is of rank one. Hence, for any T E [0,T C ), the 
minimizer of J~t in is also a minimizer of + T (3(M). From j!4|, we know that it is 
unique and given by po, in which case iu,T = iht,o — TS[po] = im,o + T f3[M]. □ 

As an immediate consequence of Lemmata 15.11 and 15.21 we obtain the following corollary. 

Corollary 5.3. Assume that (j31)-(j33) hold. There is a pure state minimizer of mass M if 
and only ifTE [0,T C ]. 

Proof. A pure state satisfies %m,t = iM,o + T (3(M) and from the concavity property stated 
in Lemma 15. II we conclude im,t < ^M,o + T (3{M) for all T > T c . □ 

We finally give a characterization of T c . 

Proposition 5.4. Assume that (/31)-(/33) hold. For any M > 0, the critical temperature 
satisfies 

T (M) - ^ ^ 

where pg and p\ are the two lowest eigenvalues of Hq defined in (|5.1I) . 
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Proof. For T < T c , there exists a unique pure state minimizer po- For such a pure state, 
the Lagrange multiplier associated to the mass constraint trpo = M is given by p = p{T). 
According to 12.71 it is given by T /3'(M) + — p{T) = for any T < T c (as long as the 
minimizer is of rank one). This uniquely determines p(T). On the other hand we know 
that / Ai = ((p? - p(T))/T) if T > (p? - /ig)//3'(M), thus proving that T c < 

(p?-p°)//?'(M). 

It remains to prove equality: By using Lemmas 15.11 and 15.21 we know that im,t c = ^M,o + 
T C (3(M). Let p be a minimizer for T = T c . The two inequalities, zm,o < £f/[p] and (3(M) < 
tr/3(p) hold as equalities if and only if, in both cases, p is of rank one. Consider a sequence 
(T( n ))„ eN such that Main^T^ = T c , > T c for any n G N and, if (p (n) ) nG N denotes a 
sequence of associated minimizers with ) £ n an< ^ < as in the proof of Lemma 15.21 
we have p^ > p^ so that > for any n G N. The sequence (p( n )) ng ^ is minimizing for 
*m,t c) thus proving that linij^oo = 0, so that lim n ^ 00 /i( n ) = p^- Passing to the limit in 

Mp^ = £ Aj n) (p$ n) + rW /3'(Aj n) ) 

completes the proof. 



□ 



6. Estimates on the maximal temperature 



All above results require T < T* , the maximal temperature. In some situations, we can 
prove that T* is finite. 

Proposition 6.1. Let /3(s) = s p with p G (1,7/5). Then, for any M > 0, the maximal 
temperature T* = T*(M) is finite. 

Proof. Let V be a given non-negative potential. From [7], we know that 

2T tr/3(p) + tr(-Ap) - tr(V» > -(2T)~F^ ( p _ ^ \pj{V)[ 1 



where 7 = and Pj(V) denotes the negative eigenvalues of —A — V . The sum is extended 
to all such eigenvalues. By the Lieb-Thirring inequality, we have the estimate 

with (7 = 7+1. In summary, this amounts to 

2T tr/3(p) + tr(-Ap) - tr(Vp) > -{2T)~^ (p - l)p~i^i C LT (7) / |^| 9 dx. 



Applying the above inequality to V = V p = n p * 



we find that 



Fr\p\ = \ tr(- A p) + I [(2 T) tr /3(p) + tr(- A p) - tr(V p p) 



> - tr(-Ap) - r p-i (2p) ^ C7 LT (7) 
Next, we invoke the Hardy-Littlewood-Sobolev inequality 

|^| 9 dx<C H Lslhp|ll. 



I^I'ds . 



/ 
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for some r > 1 such that ~ = | + ~. Notice that r > 1 means q > 3 and hence p < 3. Holder's 
inequality allows to estimate the right hand side by 

.. || || ||0 || ||-j^ Q 

\\ n p\\L r (m 3 ) < ll n pllz,i(R3) Il n plli3 



with 6 = | (i — i). Since ||n.p||^3(K3) is controlled by HVy^H^ using Sobolev's embedding, 
which is itself bounded by tr(— A p), we conclude that 

/ \V p \ q dx < cM qe {tr(-Ap)) q{1 - 9) 
for some positive constant c and, as a consequence, 
(6.1) F T [p] > \ tr(-Ap) - Ktr(-A / 9)' ? ( 1 - e \ 

for some K > 0. Moreover we find that 

7 — hp 



q (1 — 6) = 1 + rj with rj 



4(p-l) 



so that r\ is positive if p £ (1, 7/5). 

Assume that %m,t < and consider an admissible p S $)m such that J^tIp] = ^m,t- Since 
tr f3(p) is positive, as in the proof of (|2.ip . we know that for some positive constant C, which 
is independent of T > 0, 

0> F T [p] >£h[p] >tr(-Ap)-CM 3 / 2 tr(-Ap)i, 
and, as a consequence, 

tr(-A p) <C 2 M 3 . 
On the other hand, by (|6.1|) . we know that J^tIp] < means that 



tr(-Ap) > 



i 

' rr — \ >? 
1 P- 1 \ 



2K J 

The compatibility of these two conditions amounts to 

T^i < 2KC 2,? M 3 \ 

which provides an upper bound for T*{M). □ 

Finally, we infer the following asymptotic property for the infimum of J^tIp]- 

Lemma 6.2. Assume that (f31)-(f32) hold. IfT* < +oo, then limy-^* %mt = 0- 

Proof. The proof follows from the concavity of T \— )■ im,t (see Lemma 15. ip . Let pt denote 
the minimizer at To < T* , with J~t [pt ] = ~°~ f° r some 5 > 0. Then we observe 

im,t < (T - T ) ^ /3(Aj) + T To [pt ] <(T- T ) /3(M) - 5 < , 
jeN 

for all T such that: T — To < 5/(3(M), which is in contradiction with the definition of T* 
given in (|1.8p if lim inf t-> t* ^m,t < 0. □ 
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7. Concluding remarks 

Assumption (/33) is needed for Corollary 12.61 which is used itself in the proof of Proposi- 
tion (compactness of minimizing sequences). When /3(s) = s p , this means that we have 
to introduce the restriction p < 3. If look at the details of the proof, what is really needed is 
that \i = -g^jp- takes negative values. To further clarify the role of the threshold p = 3, we 
can state the following result. 

Proposition 7.1. Assume that (3(s) = s p for some p > 1. Then we have 

(7.1) M-^- + (3-p)T-^<3i M ,T 

and, as a consequence: 

(i) ifp<3, then i M ,T < (|^) 3 Uf,T /or any M > M > and T > 0. 

(ii) i/p > 3, toen i M ,T < (^) 3/(3_p) U/,T /or any M > and T > T > 0. 

Proof. Let p & Sjm and, using the representation (jl.ip . define 

With M\p] := tr p = J R3 n p dx, we find that 

M[p x ] =XM[p] = XM 

and 

As a consequence, we have 

i\M,\3-P T < A 3 !m,T , 

which proves (|7.ip by differentiating at A = 1. In case (i), since T h-> «m,t is non-decreasing, 
we have 

'AM ,T < i\M ,\ 3 -PT < A 3 «M ,T V A > 1 

and the conclusion holds with A = M/Mq. In case (ii), since M i— > im,t is non-increasing, we 
have 

«M,A3-PT < hM,\3-PT < A 3 U/,T V A G (0, 1) 

and the conclusion holds with A = (T/To) 1 / ( - 3_p ). □ 

Assume that /3(s) = for any s£l + . We observe that for T < T*(M), < jji M ,T 

if P < 3, but we have no such estimate if p > 3. In Proposition 12.31 (iii), the sufficient 
condition for showing that T*(M) = oo is precisely p > 3. Hence, at this stage, we do not 
have an example of a function (3 satisfying Assumptions (/31) and (/?2) for which existence of 
a minimizer of im,t m ?)M is granted for any M > and any T > 0. In other words, with 
T* can be infinite for a well chosen function /3, for instance /3(s) = s p , s G M + , for p > 3. 
However, in such a case we do not know if the Lagrange multiplier /i(T) is negative for any 
T > and as a consequence, the existence of a minimizer corresponding to i m,t is an open 
question for large values of T. 
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Appendix A. Proof of Proposition 13.11 

Consider the minimizer p of Proposition 13.11 and let p < be the Lagrange multiplier 
corresponding to the mass constraint tr p = M. Define 

Grip] := F t[p\ -li tr(p) . 

The density operator p is a minimizer of the unconstrained minimization problem infpgy-, Q^[p}. 
By the same argument as in the proof of Proposition 12.41 we know that p also solves the 
linearized minimization problem inf a( z^ [a] where 

£>] :=tT[{H p -p + T(3'(p))a] . 

Consider the cut-off functions XR an d £r defined in (|2.2|) and let p R := XRPXR- By 
Lemma 12.21 we know that, as R — > oo, 

tr(-Ap) >tr(-Apfl)+tr(-A(&p&i)) - 
for some positive constant C. Next we rewrite the potential energy as 

n P (x)x R (y)n p (y) , , /V X R/A {x)n p {x)^ R {y)n p {y) 



£pot[p]= ff Mx)x R (y)n P (y) dxdy+ ff 

JJm?xRZ F - 2/| JJrs 



|x - y| 



dx dy 



£|/4 (a?) n p (x) £|(y) n p (y) 



+ // ^ /4 ' ' ; - 7"' ' - dxdy . 



f - y\ 

In the second integral we use the fact that |x — y\ > R/2, whereas the third integral can be 
estimated by Lemma 13.21 Using the fact that 

e(R) := tr(-A 

= J2 X J N(tR^)\ 2 dx<2^\\Val~ + 2Y, X 4 e||V^| 2 dx 

converges to as R — > oo, we obtain that ||C^/ 4 n p * I ■ < C \/e(-R/4) — )• and can 

estimate the third integral by 

J Jr^xm. 3 \ x -y\ Jr3 

In summary this yields 

£ pot [p] < tr(y p p R ) + o(l) / £| «p dx . 
Collecting all estimates, we have proved that 



C»[pr] < m P ] - e{R) + (p + o(l)) / e R n P dx + 



C_ 

R? 



as R — > oo. Recall that e(R) is non-negative, p is negative (by Lemma 12. 5p and p is a 
minimizer of £ M so that £ M [p] < C^[p R }. As a consequence, 

(p + o(l)) / fi R n p dx + >0 
7R3 

for i? large enough, which completes the proof of Proposition 13.11 □ 
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